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Abstract 

We present new measurements of CP-violation parameters in B° — > 4>K°, K + K~Kg, fo(980)Kg, 
n'K° s , ojK° s , K%tt°, and iC° 7 (K*° -► K%ir°) decays based on a sample of 275 x 10 6 BB pairs 
collected at the T(4S I ) resonance with the Belle detector at the KEKB energy-asymmetric e + e~ 
collider. One neutral B meson is fully reconstructed in one of the specified decay channels, and the 
flavor of the accompanying B meson is identified from its decay products. CP- violation parameters 
for each of the decay modes are obtained from the asymmetries in the distributions of the proper- 
time intervals between the two B decays. All results are preliminary. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 
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I. INTRODUCTION 



The phenomena of CP violation in the flavor- changing b — > s transition are sensitive to 
physics at a very high-energy scale [1]. Theoretical studies indicate that large deviations 
from standard model (SM) expectations are allowed for time-dependent CP asymmetries in 
B° meson decays [2]. Experimental investigations have recently been launched at the two B 
factories, each of which has produced more than 10 8 BB pairs. Belle's previous measurement 
of the B° — > (pKg decay [3], which is dominated by the b — > sss transition, yielded a value 
that differs from the SM expectation by 3.5 standard deviations [4]. Measurements with a 
larger data sample are required to elucidate this difference. It is also essential to examine 
additional modes that may be sensitive to the same b — > s penguin amplitude. In this spirit, 
experimental results using the decay modes B° — > (pK®, K + K^K^, f (980)Kg, rj'Kg, and 
Kgir have already been reported [4, 5] . 

In the SM, CP violation arises from an irreducible phase, the Kobayashi-Maskawa (KM) 
phase [7], in the weak-interaction quark- mixing matrix. In particular, the SM predicts CP 
asymmetries in the time-dependent rates for B° and B° decays to a common CP eigenstate 
fcp [8]. In the decay chain T(AS) — > B°B° — > f C pj 'tag, where one of the B mesons decays 
at time tcp to a final state fcp and the other decays at time t tag to a final state /tag that 
distinguishes between B° and _B°, the decay rate has a time dependence given by 



Here S and A are CP- violation parameters, r^o is the B° lifetime, Am^ is the mass difference 
between the two B° mass eigenstates, At = t C p — t tag , and the 6-flavor charge q — +1 
(— 1) when the tagging B meson is a B° (B°). To a good approximation, the SM predicts 
S = — £/ sin 20i, where = +1(— 1) corresponds to CP-even (-odd) final states, and A = 
for both b — > ccs and b — > sss transitions. Recent measurements of time-dependent CP 
asymmetries in B° — > J/ipK® and related decay modes, which are governed by the b — ■> ccs 
transition, by Belle [9, 10] and BaBar [11] already determine sin20i rather precisely; the 
present world average value is sin20i = +0.726 ± 0.037 [6]. This serves as a firm reference 
point for the SM. 

Belle's previous measurements for B° — > <pKg, K + K~Kg and rj'Kg were based on a 140 
fb _1 data sample (DS-I) containing 152 x 10 6 BB pairs. While 4>K° S and r/'Kg final states 
are CP eigenstates with £f = —1, the K + K~K° S final state is in general a mixture of both 

= +1 and —1. Excluding K + K~ pairs that are consistent with a <fi — > K + K~ decay 
from the B° — > K + K~K§ sample, we find that the K + K~K$ state is primarily £j = +1; a 
measurement of the £/ = +1 fraction with DS-I gives 1.03 ± 0.15(stat) ± 0.05(syst) [4]. In 
the following determination of S and A, we fix = +1 for this mode. 

In this report, we describe improved measurements incorporating an additional 113 fb _1 
data sample that contains 123 x 10 6 BB pairs (DS-II) for a total of 275 x 10 6 BB pairs. 
We include additional (f)Kg and rfK® subdecay modes that were not used in the previous 
analysis. We also describe new measurements of CP asymmetries for the following CP- 
eigenstate B° decay modes: B° -> §K\ and / (980)1^ for ^ = +1; B° -> uK° s and 
KgiT for £f = — 1. The decays 5° — > <pKg and §K\ are combined in this analysis by 
redefining S as —^fS to take the opposite CP parities into account, and are collectively 
called "-B — > <pK°" . The CP asymmetries for the decay B° — > uiKg are measured for the 
first time. 



V(At) 
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1 + q ■ S sin( Am d At) + A cos( Am d At) 
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Finally, we also measure time-dependent CP violation in the decay B° — > K* 0/ ~f (K*° — > 
K s tt°) [12], which is not a CP eigenstate but is sensitive to physics beyond the SM [13]. 
Within the SM, the photon emitted from a B° (B°) meson is dominantly right-handed 
(left-handed). Therefore the polarization of the photon carries information on the origi- 
nal 6-flavor; the decay is thus almost flavor-specific. The SM predicts a small asymmetry 
S ~ — 2(m s /mfo)sin2(/> 1 , where m b {m s ) is the 6-quark (s-quark) mass [13]. Any significant 
deviation from this expectation would be a manifestation of physics beyond the SM. 

At the KEKB energy- asymmetric e + e~ (3.5 on 8.0 GeV) collider [14], the T(4S') is pro- 
duced with a Lorentz boost of (3^ = 0.425 nearly along the electron beamline (z). Since 
the B° and B° mesons are approximately at rest in the T(4S) center-of-mass system (cms), 
At can be determined from the displacement in z between the fcp and /tag decay vertices: 
At ~ {z CP - z tag )/{^c) = Az/(l3>yc). 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside of the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [15]. Two inner detector configurations 
were used. A 2.0 cm radius beampipe and a 3-layer silicon vertex detector (SVD-I) were 
used for DS-I, while a 1.5 cm radius beampipe, a 4-layer silicon detector (SVD-II) and a 
small-cell inner drift chamber were used for DS-II [16]. 

II. EVENT SELECTION, FLAVOR TAGGING AND VERTEX RECONSTRUC- 
TION 

A. Overview 

We reconstruct the following B° decay modes to measure CP asymmetries: B° — > 4>K S , 
<pK° L , K + K-K° s , f (980)K s , r)'K s , uK° s , K s n°, and K*°-f (K*° -> K s n°). We exclude 
K + K~ pairs that are consistent with a <fi — > K + K~ decay from the B° — > K + K~K S sample. 
The intermediate meson states are reconstructed from the following decays: n° — > 77, 
Kg — > 7r + 7r~ (also 7r°7r° for the (f>K s decay), r\ — > 77 or 7r + 7r~7r°, p° — > 7r + 7r~ , uj — > 7r + 7r _ 7r , 
K*° -> K S 7T°, rj' -> p°7 or ^tt+tt-, / (980) -> tt+tt-, and -> 

Among the decay chains listed above, B° — > (ifg — > 7r + 7r _ ), — > K + K~K S , 

B° — > ^'-f^g (V — > P°t), and 5° — >■ 77^5 I 7 ?' ~~ *■ i] 7T+7l ~i V ~ * 77) decays were used in 
the previous analysis [4]. The selection criteria for these decays remain the same. For the 
other B° decay modes that are now included, identification of photons, neutral and charged 
kaons, and neutral and charged pions is based on the procedure used previously. However, 
the selection criteria for each B° decay mode were optimized individually and are thus 
different from one another. 

B. B° -> cj>K% and K+K~K% 

We use well-reconstructed charged tracks with a sufficient number of associated hits in 
the CDC. Charged kaons and pions except for those from K s — > 7r + 7r~ decays are required 
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to originate from the interaction point (IP). We distinguish charged kaons from pions based 
on a kaon (pion) likelihood £k(it) derived from the TOF, ACC and dE/dx measurements in 
the CDC. 

Pairs of oppositely charged tracks that have an invariant mass within 0.030 GeV/c 2 of the 
nominal Kg mass are used to reconstruct Kg — > 7r + 7r~ decays. The 7r + 7r~ vertex is required 
to be displaced from the IP by a minimum transverse distance of 0.22 cm for high momentum 
(> 1.5 GeV/c) candidates and 0.08 cm for those with momentum less than 1.5 GeV/c. The 
direction of the pion pair momentum must also agree with the direction defined by the IP 
and the vertex displacement within 0.03 rad for high-momentum candidates, and within 0.1 
rad for the remaining candidates. 

Photons are identified as isolated ECL clusters that are not matched to any charged 
track. To select Kg — > 7r 7r decays, we reconstruct ir° candidates from pairs of photons with 
E 1 > 0.05 GeV, where E^ is the photon energy measured with the ECL. The reconstructed 
7T° candidate is required to have an invariant mass between 0.08 and 0.15 GeV/c 2 and a 
momentum above 0.1 GeV/c. The large mass range is used to achieve a high reconstruction 
efficiency Candidate Kg — > 7r°7r° decays are required to have invariant masses between 0.47 
GeV/c 2 and 0.52 GeV/c 2 , where we perform a fit with constraints on the K s vertex and the 
7T° masses to improve the 7r°7r° invariant mass resolution. We also require that the distance 
between the IP and the reconstructed Kg decay vertex be larger than —10 cm, where the 
positive direction is defined by the Kg momentum. 

Candidate — > K + K~ decays are required to have a K + K~ invariant mass that is 
within 0.01 GeV/c 2 of the nominal meson mass. Since the meson selection is effective 
in reducing background events, we impose only minimal kaon-identification requirements; 
TZk/tt = £-k/{£-k + £tt) > 0.1 is required, where the kaon likelihood ratio TZr/tt has values 
between (likely to be a pion) and 1 (likely to be a kaon). We use a more stringent kaon- 
identification requirement, IZx/n > 0.6, to select non-resonant K + K~ candidates for the 
decay B° -> K + K~K° S . We reject K + K~ pairs that are consistent with D° -> K + K~, 
Xco —> K + K~~, or J/ip — > K + K~ decays. We also remove D + — > KgK + candidates. 

For reconstructed B — > fcp candidates, we identify B meson decays using the en- 
ergy difference AE = E B ms — E^ m and the beam-energy constrained mass M^ c = 

\J(.Eblfm) 2 ~ (p c b 1s ) 2 i where E^ m is the beam energy in the cms, and E B ms and p c ^ s are 
the cms energy and momentum of the reconstructed B candidate, respectively. The B 
meson signal region is defined as \AE\ < 0.06 GeV for B° — > <pK s (K s — > 7r + 7r~), 
-0.15 GeV < AE < 0.1 GeV for B° -> <pK° s (K s -> 7r°7r°), \AE\ < 0.04 GeV for 
B° -> K + K~Kg, and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 for all decays. 

The dominant background to the B° — > (f)Kg decay comes from e + e~ — > uu, dd, ss, 
or cc continuum events. Since these tend to be jet-like, while the signal events tend to 
be spherical, we use a set of variables that characterize the event topology to distinguish 
between the two. We combine S±, $t and modified Fox- Wolfram moments [17] into a Fisher 
discriminant J 7 , where S± is the scalar sum of the transverse momenta of all particles outside 
a 45° cone around the candidate meson direction divided by the scalar sum of their total 
momenta, and 9t is the angle between the thrust axis of the B candidate and that of the 
other particles in the cms. We also use the angle of the reconstructed B° candidate with 
respect to the beam direction in the cms (9b), and the helicity angle 9h defined as the 
angle between the B° meson momentum and the daughter K + momentum in the meson 
rest frame. We combine cos 9b and cos9h into a signal [background] likelihood variable, 
which is defined as £ S i g [bk g ] = £ S i g [bk g ] (.T 7 ) x £ si g[bkg](cos0 B ) x £ s ig[bkg](cos0#). We impose 
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requirements on the likelihood ratio 7Z s /b = £ S ig/ (£si g + £bk g ) to maximize the figure-of- 

merit (FoM) defined as N^ g c /^N^ g c + N hkg , where A^ g c (N hkg ) represents the expected 

number of signal (background) events in the signal region. We estimate N^ g using Monte 
Carlo (MC) events, while N hkg is determined from events outside the signal region. The 
requirement for H s /b depends both on the decay mode and on the flavor-tagging quality, r, 
which is described in Sec. III. The threshold values range from 0.1 (used for r > 0.875) to 
0.4 (used for r < 0.25) for the decay B° -> <j>K° s (K° s -> tt+tt-), and from 0.25 to 0.65 for 
the decay B° — > K + K~K S . We impose a more stringent requirement, TZ s /b > 0.75, for all r 
values in the decay B° -> <pK° s (K° s -> tt tt°). 

We use events outside the signal region as well as a large MC sample to study the 
background components. The dominant background is from continuum. The contributions 
from BB events are small. The contamination of B° — > K + K~K S events in the B° — > <pK s 
sample is 7.1 ± 1.6% (6.2 ± 2.0%) for DS-I (DS-II). Backgrounds from the decay B° -> 
fo(980)K s (/ (980) — > K + K~), which has a CP eigenvalue opposite to (f)K s , are found to 
be 0.4i^% (0.0±g:g%) for DS-I (DS-II). The influence of these backgrounds is treated as a 
source of systematic uncertainty. 

Figure 1(a) and (c) show the M bc distribution for the reconstructed B° — > <pK s and 
K + K~K S candidates within the AE signal regions after flavor tagging and vertex recon- 
struction. The signal yield is determined from an unbinned two-dimensional maximum- 
likelihood fit to the AE-Mh c distribution in the fit region defined as 5.2 GeV/c 2 < 
M bc < 5.3 GeV/c 2 for all modes, and -0.12 GeV < AE < 0.25 GeV for the B° -> 
(f)K° s (K° s -> tt+tt-) or K + K-K° s decay and -0.15 GeV < AE < 0.25 GeV for the 
B° -> (pK° s (K° s -> 7r 7T°) decay. The <j>K° s (K° s -> tt+tt-) signal distribution is mod- 
eled with a Gaussian function (a sum of two Gaussian functions) for M\, c (AE). The 
4>Kg (Kg — > tt tt°) signal distribution is modeled with a smoothed histogram obtained from 
MC events. For the continuum background, we use the ARGUS parameterization [18] for 
M bc and a linear function for AE. The fits yield 139 ± 14 B° -> <pK s events and 399 ± 28 
B° — > K + K~K S events in the signal region, where the errors are statistical only. 

C. B° -> cj)Kl 

Candidate <fi — > K + K~ decays are selected with the same criteria as described above. 
We select K\ candidates based on KLM and ECL information. There are two classes of 
K\ candidates that we refer to as KLM and ECL candidates. The requirements for the 
KLM candidates are the same as those used in the B° — > J/ipK^ selection for the sin20! 
measurement [9]. ECL candidates are selected from ECL clusters using a K\ likelihood 
ratio [9], which is calculated from the following information: the distance between the ECL 
cluster and the closest extrapolated charged track position; the ECL cluster energy; E 9 /E 2 5, 
the ratio of energies summed in 3 x 3 and 5x5 arrays of CsI(Tl) crystals surrounding the 
crystal at the center of the shower; the ECL shower width and the invariant mass of the 
shower. The likelihood ratio is required to be greater than 0.8. For both KLM and ECL 
candidates, we also require that the cosine of the angle between the K\ direction and the 
direction of the missing momentum of the event in the laboratory frame be greater than 0.6. 

Since the energy of the K\ is not measured, M bc and AE cannot be calculated in the 
same way that is used for the other final states. Using the four-momentum of a reconstructed 
candidate and the K\ flight direction, we calculate the momentum of the K\ candidate 
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requiring AE = 0. We then calculate p§" s , the momentum of the B candidate in the cms, 
and define the B meson signal region as 0.2 GeV/c < p c ™ < 0.5 GeV/c. We impose the 
requirement 1Z s /b > 0.98 to reduce the continuum background. Here 7Z s /b is based on the 
discriminating variables used for the B° — > <pK s decay and the number of tracks originating 
from the IP with a momentum above 0.1 GeV/c. The TZ s /b requirement is chosen to optimize 
the FoM, which is calculated taking the background from both continuum and generic B 
decays into account. The K\ detection efficiency difference between data and MC is studied 
using the decay B° — > J/ipK®, and corrections are applied to the B° — > MC events 
to calculate the FoM. If there is more than one candidate B° — > 4>K^ decay in the signal 
region, we take the one with the highest 1Z s /b value. ECL candidates are not used if there 
is a candidate B° — > (f)Kl decay with a KLM candidate. We find that about 90% of signal 
events are reconstructed with KLM candidates. 

We study the background components using a large MC sample as well as data taken 
with cms energy 60 MeV below the nominal T(4S) mass (off-resonance data). The dominant 
background is from continuum. A MC study with the efficiency correction obtained from 
B° — > J/ipK^ data yields 9±5 background events from B decays, which include B° — > (j)K*°, 
4>K S and B + — > <pK* + decays. The influence of these backgrounds, including their CP 
asymmetries, is treated as a source of systematic uncertainty. 

The p c g ns distribution after flavor tagging and vertex reconstruction is shown in Fig. 1(b). 
The signal yield is determined from an extended unbinned maximum-likelihood fit in the 
range GeV/c < p c ™ < 1 GeV/c. The B° — > §K\ signal shape is obtained from MC 
events. Background from BB pairs is also modeled with MC. We fix the ratio between the 
signal yield and the BB background based on known branching fractions and reconstruction 
efficiencies; uncertainty in the ratio is treated as a source of systematic error. The continuum 
background distribution is represented by a smoothed histogram obtained from MC events; 
we confirm that the function describes the off-resonance data well. The fit yields 36 ±15 ±10 
B° — > 4>Kl events, where the first error is statistical and the second error is systematic. The 
sources of the systematic error include uncertainties in the efficiency corrections, in BB 
background branching fractions and in the background parameterizations. The result is in 
good agreement with the expected B° — > (fiK® signal yield (36 ± 9 events) obtained from 
MC after applying the efficiency correction from the B° — > J/i/jK® data. 

D. B° -»• / (980)K° 

Candidate K s — > 7r + 7r~ decays are selected with the criteria that are slightly different from 
those used for the B° — > <pK s decay to obtain the best performance for the B° — > f (980)K s 
decay. Pairs of oppositely charged pions that have invariant masses between 0.890 and 1.088 
GeV/c 2 are used to reconstruct / (980) — > it + it~ decays. Tracks that are identified as kaons 
(TZr/tt > 0.7) or electrons are not used. We require that both K s n + and K s n~ combinations 
have invariant masses more than 0.1 GeV/c 2 above the nominal charged D meson mass; this 
removes background from — > K^tt^ and K*^ — * decays. 

The B meson signal region is defined as \AE\ < 0.06 GeV and 5.27 GeV/c 2 < M bc < 
5.29 GeV/c 2 . The dominant background is from continuum. For the continuum suppression, 
we require 7Z s /h > 0.6 for events with the best-quality flavor tagging (r > 0.875), and 
7Z s /b > 0.8 for other events. Here the signal likelihood ratio 7Z s /b is obtained from cos 9b 
and J 7 , which consists of the modified Fox- Wolfram moments and cos 9 T . 

Figure 1(d) shows the M bc distribution for the reconstructed B° — > f (980)K s candidates 
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within the AE signal region after flavor tagging and vertex reconstruction. For the signal 
yield extraction, we first perform an unbinned two-dimensional maximum-likelihood fit to 
the AE-M hc distribution in the fit region defined as 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 
and —0.3 GeV < AE < 0.4 GeV. The signal is modeled with a Gaussian function (a 
sum of two Gaussian functions) for M^ c {AE). For the continuum background, we use the 
ARGUS parameterization for M^ c and a linear function for AE. The fit yields the number 
of B° — > ty + ti^K s events that have invariant masses within the / (980) resonance 

region, which may include contributions from B° — > p°K s as well as non-resonant three-body 
B° — > n + 7r~Kg decays. To separate these peaking backgrounds from the B° — > f {980)K s 
decay, we perform another fit to the 7i + ir~ invariant mass distribution for the events inside 
the AE-M^c signal region. We use Breit-Wigner functions for the B° — > fo{980)K s signal as 
well as for B° — > pK s and a possible resonance above the / (980) mass region [19]. Three- 
body B° — > it + it~K s decays are modeled with a fourth-order polynomial function. The other 
background is modeled with a threshold function. The fit yields 94 ± 14 B° — > fo{980)K s 
events. The peaking background contribution in the AE-M^ C signal region is estimated to 
be 9 ± 3 events. 

E. B° -> <q'K% 

Candidate K s — > it + it~ decays are selected with the same criteria as those used for 
the B° — > <pK s decay. Charged pions from the rj, p° or 7/ decay are selected from tracks 
originating from the IP. We reject kaon candidates by requiring TZk/tt < 0.9. Candidate 
photons from ir° — > 77 decays are required to have E 1 > 0.05 GeV. The reconstructed 7r° 
candidate is required to satisfy 0.118 GeV/c 2 < M 77 < 0.15 GeV/c 2 and p c J s > 0.1 GeV/c, 
where M 77 and p™ s are the invariant mass and the momentum in the cms, respectively. 
Candidate photons from 77 — > 77 (7/ p 0/ ~f) decays are required to have E 1 > 0.05 (0.1) GeV. 
The invariant mass of the photon pair is required to be between 0.5 and 0.57 GeV/c 2 for the 
V ~^ 77 decay. The 7r + 7r~7r° invariant mass is required to be between 0.535 and 0.558 GeV/c 2 
for the 77 — > 7T + 7r~7r decay. A kinematic fit with an 77 mass constraint is performed using 
the fitted vertex of the 7r + 7r~ tracks from the 7/ as the decay point. For 7/ — > p 0/ y decays, 
candidate p° mesons are reconstructed from pairs of vertex-constrained 7r + 7r~ tracks with 
an invariant mass between 0.55 and 0.92 GeV/c 2 . The rj' — > rin + n^ candidates are required 
to have a reconstructed mass between 0.94 and 0.97 GeV/c 2 (0.95 and 0.966 GeV/c 2 ) for 
the rj — > 77 {rj — > 7r + 7r _ 7r°) decay. Candidate 77' — > p°7 decays are required to have a 
reconstructed mass from 0.935 to 0.975 GeV/c 2 . 

The B meson signal region is defined as \ AE\ < 0.06 GeV for B° — > rj'K s {rj' — > p°7), —0.1 
GeV < AE < 0.08 GeV for B° -> 77'^ (77' -> t^tt+tt-, 77 -> 77) or -0.08 GeV < AE < 0.06 
GeV for B° -> (V ^ ^vt+tt-, 77 -> tt+tt-tt ), and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 

for all decays. The continuum suppression is based on the likelihood ratio 7Z s /b obtained from 
the same discriminating variables used for the B° — > <pK s decay, except that cos 6h for B° 
decays is not used; we only use the helicity angle for the decay rj' — > p7 (p — > 7r + 7r _ ), which 
is defined as the angle between the rj' meson momentum and the daughter ir + momentum in 
the p meson rest frame. The minimum 7Z s /h requirement depends both on the decay mode 
and on the flavor-tagging quality, and ranges from (i.e. no requirement) to 0.4. 

We use events outside the signal region as well as a large MC sample to study the back- 
ground components. The dominant background is from continuum. In addition, according 
to MC simulation, there is a small (~ 3%) contamination from BB background events in 
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B° — > 77^5 (77' — > p°7)- The contributions from BB events are smaller for other modes. 
The influence of these backgrounds is treated as a source of systematic uncertainty. 

Figure 1(e) shows the M bc distribution for the reconstructed B° — > r)'K s candidates 
within the AE values signal region after flavor tagging and vertex reconstruction. The 
signal yield is determined from unbinned two-dimensional maximum-likelihood fits to the 
AE-M\j C distributions in the fit region defined as 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and 
—0.25 GeV < AE < 0.25 GeV. We perform the fit for each final state separately. The r)'K s 
signal distribution is modeled with a sum of two (three) Gaussian functions for M bc (AE). 
For the continuum background, we use the ARGUS parameterization for M bc and a linear 
function for AE. For the rj' — > p7 mode, we include the BB background shape obtained 
from MC in the fits. The fits yield a total of 512 ±27 B° — > r]'K s events in the signal region, 
where the error is statistical only. 

F. B° -> uK% 

Candidate Kg — > it + it~ decays are selected with criteria that are identical to those used 
for the B° — > <pK s decay. Pions for the uo — > 7r + 7r~7r° decay are selected with the same 
criteria used for the 77 — > 7r + 7r~7r° decay, except that we require p"o s > 0.35 GeV/c. The 
7T + 7r~7r° invariant mass is required to be within 0.03 GeV/c 2 of the nominal uo mass. The B 
meson signal region is defined as \AE\ < 0.06 GeV and 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 . 
The dominant background is from continuum. The continuum suppression is based on the 
likelihood ratio 7Z s /b obtained from the same discriminating variables used for the B° — > <f)K s 
decay; the helicity angle Qu is defined as the angle between the B° meson momentum and 
the cross product of the 7i + and n~ momenta in the uj meson rest frame. The minimum 
7l s /h requirement depends both on the decay mode and on the flavor-tagging quality, and 
ranges from 0.3 (used for r > 0.875) to 0.9 (used for r < 0.25). The contribution from BB 
events is negligibly small. 

Figure 1(f) shows the M bc distribution for the reconstructed B° — > uK s candidates 
within the AE signal region after flavor tagging and vertex reconstruction. The sig- 
nal yield is determined from an unbinned two-dimensional maximum-likelihood fit to the 
AE-M^c distribution in the fit region defined as 5.2 GeV/c 2 < M bc < 5.3 GeV/c 2 and 
—0.12 GeV < AE < 0.25 GeV. The AE distribution is modeled with a sum of two (three) 
Gaussian functions for M bc (AE). For the continuum background, we use the ARGUS 
parameterization for M bc and a linear function for AE. The fit yields 31 ± 7 B° — > ujK s 
events in the signal region with a statistical significance (X) of 7.3, where X is defined as 

E = ^—2 \n(C /CN sig ), and C and ^7v sig denote the maximum likelihoods of the fits without 
and with the signal component, respectively. 

G. B° -> K° S 7V° 

Candidate K s — > it + it~ decays are selected with the same criteria as those used for the 
B° — > <$K S decay, except that we impose a more stringent invariant mass requirement; only 
pairs of oppositely charged pions that have an invariant mass within 0.015 GeV/c 2 of the 
nominal K s mass are used. The n° selection criteria are the same as those used for the 
B° -> rfK° s decay. 
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The B meson signal region is defined as —0.15 GeV < AE < 0.1 GeV and 5.27 GeV/c 2 < 
Mbc < 5.29 GeV/c 2 . The dominant background is from continuum. We use extended 
modified Fox- Wolfram moments, which were applied for the selection of the B° — > 7r°7r° 
decay [20], to form F . We then combine likelihoods for T and cos9 B to obtain the event 
likelihood ratio TZ s /h for continuum suppression. As described below, we include events that 
do not have B decay vertex information in our fit to obtain a better sensitivity for the CP- 
violation parameter A. For events with vertex information, the high-7?. s /b region is defined 
as TZ s/h > 0.74 (0.76) for DS-I (DS-II), and the low-ft s/b region as 0.4 < TZ s/h < 0.74 (0.76) 
for DS-I (DS-II). For events without vertex information, the high-7?. s /b region is defined as 
U s/h > 0.78 and the low-TZ s/b region as 0.4 < TZ s/h < 0.78 for both DS-I and DS-II. 

Figure 1(g) shows the M bc distribution for the high-7?. s /b B° — > K s tt° candidates within 
the AE signal region after flavor tagging and before vertex reconstruction. Also shown in 
Fig. 1(h) is the M bc distribution for the low-7?. s /b B° — > K s n° candidates. The signal yield 
is determined from an unbinned two-dimensional maximum-likelihood fit to the AE-M bc 
distribution in the fit region defined as 5.2 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.2 GeV < 
AE < 0.5 GeV. The B° — > K s n° signal distribution is modeled with a Gaussian function 
for M bc and with a Crystal Ball function for AE. For the continuum background, we use 
the ARGUS parameterization for M bc and a second-order Chebyshev function for AE. The 
BB background is negligibly small and its influence is treated as a source of systematic 
uncertainty. The fits yield 168 ± 16 and 83 ± 18 B° — > K s n events in the high-7?. s /b 
and low-7£ s /b signal regions, respectively, where the errors are statistical only. The same 
procedure after the vertex reconstruction yields a total of 77 ± 13 K s ir° events. 

H. B° -> K*° 1 (K*° -> K°7T°) 

The selection criteria are optimized for the measurement of time-dependent CP asym- 
metries and are thus different from those used in Belle's previous measurement of B — > K*^ 
branching fractions and direct CP asymmetries [21]. Candidate K s — > it + it~ decays are 
selected with the same criteria as those used for the B° — > 4>K S decay, except that we im- 
pose a more stringent invariant mass requirement; only pairs of oppositely charged pions 
that have an invariant mass within 0.006 GeV/c 2 of the nominal K s mass are used. The 
7T° selection criteria are the same as those used for the B° — » r)'K s decay except for a more 
stringent ir° momentum requirement, p c ™ > 0.3 GeV/c. The K s ir° invariant mass, M K o n o, 
is required to be between 0.6 and 1.8 GeV/c 2 . 

Prompt photons from the -B — > K* 0/ ~f decay are required to satisfy 1.4 GeV < E^ ms < 
3.4 GeV, where E^ ms is the photon energy in the cms. If there is more than one candidate, 
the one with the largest E^ ms is selected. For the selected photon, we require E 9 /E 2 5 > 0.95, 
where E 9 /E 2 5 is defined in Sec. II C. Photons for candidate tc° — > 77 or rj — >■ 77 decays are 
not used; we reject photon pairs that satisfy £ n o > 0.18 or C v > 0.18, where C n o^ is a n° 
(77) likelihood described in detail elsewhere [22]. The polar angle of the photon direction in 
the laboratory frame is required to be between 33° and 128° for DS-I, while no requirement 
is imposed for DS-II as the material within the acceptance of the ECL is much reduced for 
this dataset. 

Candidate K* + — > K s 7r + decays are also selected using a similar procedure to reconstruct 
the decay B° -> K* + -f. Candidate B° -> K*°-f (K*° -> K s ir°) and B + -> K* + ^ (K* + -> 
K s n + ) decays are selected simultaneously; we allow only one candidate for each event. The 
best candidate selection is based on the event likelihood ratio 1Z s / h that is obtained by 
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combining J 7 , which uses the extended modified Fox- Wolfram moments as discriminating 
variables, with cos 8h defined as the angle between the B° meson momentum and the daugh- 
ter Kg momentum in the K*° meson rest frame. We select the candidate with the largest 

The signal region for the B° — > K*°*y (K*° — > K s n°) decay is defined as —0.2 GeV < 
AE < 0.1 GeV, 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and 0.8GeV/c 2 < M K o n o < 1.0 GeV/c 2 . 
We require 7Z s /b > 0.5 to reduce the continuum background. 

We use events outside the signal region as well as a large MC sample to study the back- 
ground components. The dominant background is from continuum. Background contribu- 
tions from B decays are significantly smaller than those from continuum, and are dominated 
by cross-feed from B + — > K* +/ ~f decays, other radiative B decays and charmless B decays. 
Background from other BB decays is found to be negligible. 

Figure l(i) shows the M bc distribution for the reconstructed B° — > K* ^ candidates 
within the AE signal region after flavor tagging and vertex reconstruction. The signal yield 
is determined from an unbinned two-dimensional maximum-likelihood fit to the AE-M^ C 
distribution in the fit region defined as 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.5 GeV < 
AE < 0.5 GeV. The B° — > K* ^ signal distribution is represented by a smoothed histogram 
obtained from MC simulation that accounts for the correlation between M bc and AE. The 
background from B decays is also modeled with a smoothed histogram obtained from MC 
events; its normalization is a free parameter in the fit. For the continuum background, we 
use the ARGUS parameterization for M bc and a second-order Chebyshev function for AE. 
The fit yields 57 ± 9 B° -> K*°j (K*° -> K s ir°) events, where the error is statistical only. 
For reference, we also measure the signal before vertex reconstruction and obtain 132 ± 14 
events. 

I. Flavor Tagging 

The 6-flavor of the accompanying B meson is identified from inclusive properties of par- 
ticles that are not associated with the reconstructed B° — > f C p decay. We use the same 
procedure that is used for the sin 20! measurement [10]. The algorithm for flavor tagging is 
described in detail elsewhere [23]. We use two parameters, q and r, to represent the tagging 
information. The first, q, is already defined in Eq. (1). The parameter r is an event-by- 
event, MC-determined flavor-tagging dilution factor that ranges from r = for no flavor 
discrimination to r = 1 for unambiguous flavor assignment. It is used only to sort data into 
six r intervals. The wrong tag fractions for the six r intervals, wi (I = 1, 6), and differences 
between B° and B° decays, Awi, are determined from the data; we use the same values 
that were used for the sin 20! measurement [10] for DS-I. Wrong tag fractions for DS-II 
are separately obtained with the same procedure; we find that the values for DS-II, which 
are listed in Table I, are slightly smaller than those for DS-I. The total effective tagging 
efficiency for DS-II is determined to be e e g = Z)f=i e «(l — 2w/) 2 = 0.30 ± 0.01, where q is 
the event fraction for each r interval determined from the J/ipK s data and listed in Table I. 
The error includes both statistical and systematic uncertainties. 
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J. Vertex Reconstruction 

The vertex position for the fcp decay is reconstructed using charged tracks that have 
enough SVD hits. A constraint on the IP is also used with the selected tracks; the IP profile 
is convolved with finite B flight length in the plane perpendicular to the z axis. The pions 
from Kg decays are not used except for the analysis of B° — > K* 0/ ~f and B° — > K s tt° decays. 
The vertex for B° — > K* 0/ ~f and B° — » K s 7i° decays is reconstructed using the K s trajectory 
and the IP constraint, where both pions from the K s decay are required to have enough 
SVD hits to reconstruct a vertex. The reconstruction efficiency depends both on the K s 
momentum and on the SVD geometry; efficiencies with SVD-II are higher than those with 
SVD-I because of the larger outer radius and the additional layer. 

The /tag vertex determination with SVD-I remains unchanged from the previous pub- 
lication [4], and is described in detail elsewhere [24]; to minimize the effect of long-lived 
particles, secondary vertices from charmed hadrons and a small fraction of poorly recon- 
structed tracks, we adopt an iterative procedure in which the track that gives the largest 
contribution to the vertex x 2 is removed at each step until a good x 2 is obtained. 

For SVD-II, we find that the same vertex reconstruction algorithm results in a larger 
outlier fraction when only one track remains after the iteration procedure. Therefore, in 
this case, we repeat the iteration procedure with a more stringent requirement on the SVD- 
II hit pattern. The resulting outlier fraction is comparable to that for SVD-I, while the 
inefficiency caused by this change is small (2.5%). 

K. Summary of Signal Yields 

The signal yields for B° — > fcp decays, iV B i g , after flavor tagging and vertex reconstruction 
(before the vertex reconstruction for the decay B° — > K s tt°) are summarized in Table II. 
The signal purities are also listed in the table. 

III. RESULTS OF CP ASYMMETRY MEASUREMENTS 

We determine S and A for each mode by performing an unbinned maximum-likelihood 
fit to the observed At distribution. The probability density function (PDF) expected for 
the signal distribution, V s i g (At; S, A, q, wi, Aw{), is given by Eq. (1) incorporating the effect 
of incorrect flavor assignment. The distribution is convolved with the proper-time interval 
resolution function R sig (At), which takes into account the finite vertex resolution. 

For the decays B° -> <j>K%, K + K~K° S , <j>K° L , f (980)K s , r)'K° s and uK° s , we use flavor- 
specific B decays governed by semileptonic or hadronic b — > c transitions to determine 
the resolution function. We perform a simultaneous multi-parameter fit to these high- 
statistics control samples to obtain the resolution function parameters, wrong-tag fractions 
(Section III), Am^, r B + and r B o. We use the same resolution function used for the sin20i 
measurement for DS-I [10]. For DS-II, the following modifications are introduced: a sum 
of two Gaussian functions is used to model the resolution of the fcp vertex while a single 
Gaussian function is used for DS-I; a sum of two Gaussian functions is used to model the 
resolution of the tag-side vertex obtained with one track and the IP constraint, while a 
single Gaussian function is used for DS-I. These modifications are needed to account for 
differences for SVD-I and SVD-II, as well as different background conditions in DS-I and 
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DS-II. We test the new resolution parameterization using MC events on which we overlay 
beam-related background taken from data. A fit to the MC sample yields correct values for 
all parameters. With the multi-parameter fit to data, we find that the standard deviation of 
the main Gaussian component of the resolution function is reduced from 78 /im to 55 /xm, 
which is consistent with our expectation from the improved impact parameter resolution of 
SVD-II [16]. The same fit also yields t b o = 1.518 ± 0.012 ps, r B + = 1.652 ± 0.014 ps and 
Arrid = 0.516 ± 0.007 ps -1 . The results are consistent with those obtained with DS-I [10] 
and also with the world average values [26] . Thus we conclude that the resolution of SVD-II 
is well understood. 

For the decays B° — > Kgir and K* ^ (K*° — > Kg-rr ), we use the resolution function 
described above with additional parameters that rescale vertex errors. The rescaling pa- 
rameters depend on the detector configuration (SVD-I or SVD-II), SVD hit patterns of 
charged pions from the K s decay, and K s decay vertex position in the plane perpendicular 
to the beam axis. These parameters are determined from a fit to the At distribution of 
B° — > J/ipKg data. Here the K s and the IP constraint are used for the vertex reconstruc- 
tion, the B° lifetime is fixed at the world average value, and 6-flavor tagging information is 
not used so that the expected PDF is an exponential function convolved with the resolution 
function. 

We check the resulting resolution function by also reconstructing the vertex with leptons 
from J /if) decays and the IP constraint. We find that the distribution of the distance 
between the vertex positions obtained with the two methods is well represented by the 
obtained resolution function convolved with the well-known resolution for the J/ip vertex. 
Finally, we also perform a fit to the B° — > J/ipK s sample with 6-flavor information and 
obtain Sj/^ K o = +0.68 ± O.lO(stat) and Aj^ K o = +0.02 ± 0.04(stat), which are in good 
agreement with the world average values. Thus, we conclude that the vertex resolution for 
the B° -> K s tt° and B° -> fT*° 7 {K*° -> K s ir°) decays is well understood. 

We determine the following likelihood value for each event: 

Pi = (l-/oi) /[/ S i g Psig(AO^sig(At i -At / ) 

+ (1 - /sig)nk g (A0i4k g (A^ - At')p(At') 

+ /oiPoi(A^) (2) 

where P \(At) is a broad Gaussian function that represents an outlier component with a 
small fraction / ol . The signal probability / sig depends on the r region and is calculated on 
an event-by-event basis as a function of p ™ for the B° — > <pK L decay and as a function of 
AE and M^ c for the other modes. A PDF for background events, Pbkg(At), is modeled as a 
sum of exponential and prompt components, and is convolved with a sum of two Gaussians 
-Rbkg- All parameters in Pbk g (At) and i?bk g are determined by the fit to the At distribution of 
a background-enhanced control sample [25]; i.e. events outside of the AE-M^ C signal region. 
We fix t b o and Am^ at their world-average values [26]. In order to reduce the statistical 
error on A, we include events without vertex information in the analysis of B° — > K s 7i°. 
The likelihood value in this case is obtained by integrating Eq. (2) over Atj. 

The only free parameters in the final fit are S and A, which are determined by maximizing 
the likelihood function L = ^ P^Atf, S , A) where the product is over all events. Since A 
for the decay -B — > K* 0/ ~f is well measured in the B° — > K* 0/ ~f (K*° — > K + n~) mode and is 
consistent with zero, we fix A at zero and perform a fit to the B° — > K* 0/ ~f (K*° — * K s ir°) 
sample with S as the only free parameter. 
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Table III summarizes the fit results of S and A. We define the raw asymmetry in each 
At bin by (N q=+ \ — N q= _i) / (N q=+ i + N q= _i), where N q=+ u_u is the number of observed 
candidates with q = +1(— 1) [27]. Figures 2(a-g) show the raw asymmetries in two regions of 
the flavor-tagging parameter r. While the numbers of events in the two regions are similar, 
the effective tagging efficiency is much larger and the background dilution is smaller in the 
region 0.5 < r < 1.0. Note that these projections onto the At axis do not take into account 
event-by-event information (such as the signal fraction, the wrong tag fraction and the vertex 
resolution), which is used in the unbinned maximum-likelihood fit. 

Tables IV and V list the systematic errors on S and A, respectively. The total systematic 
errors are obtained by adding each contribution in quadrature, and are much smaller than 
the statistical errors for all modes. 

To determine the systematic error that arises from uncertainties in the vertex recon- 
struction, the track and vertex selection criteria are varied to search for possible systematic 
biases. Small biases in the Az measurement are observed in e + e~ — > fi + fi~ and other con- 
trol samples. Systematic errors are estimated by applying special correction functions to 
account for the observed biases, repeating the fit, and comparing the obtained values with 
the nominal results. The systematic error due to the IP constraint in the vertex reconstruc- 
tion is estimated by varying (±10//m) the smearing used to account for the B flight length. 
Systematic errors due to imperfect SVD alignment are determined from MC samples that 
have artificial mis-alignment effects to reproduce impact-parameter resolutions observed in 
data. 

Systematic errors due to uncertainties in the wrong tag fractions are studied by varying 
the wrong tag fraction individually for each r region. Systematic errors due to uncertainties 
in the resolution function are also estimated by varying each resolution parameter obtained 
from data (MC) by ±lcr (±2<r), repeating the fit and adding each variation in quadrature. 
Each physics parameter such as r B o and Am^ is also varied by its error. A possible fit bias 
is examined by fitting a large number of MC events. 

Systematic errors from uncertainties in the background fractions and in the background 
At shape are estimated by varying each background parameter obtained from data (MC) 
by ±lo" (±2(j). Uncertainties in the background B decay model are also considered for the 
B° — > fT*°7 (K*° — > Kg7t°) mode; we compare different theoretical models for radiative B 
decays and take the largest variation as the systematic error. 

Additional sources of systematic errors are considered for B decay backgrounds that are 
neglected in the PDF. We consider uncertainties both in their fractions and CP asymmetries; 
for modes that have non- vanishing CP asymmetries, we conservatively vary the CP- violation 
parameters within the physical region and take the largest variation as the systematic error. 
The effect of backgrounds from K + K~K° S and f (980)K° s (/ (980) -> K+K~) in the B° -> 
4>Kg sample is considered. Uncertainties from B — > <pK* and other rare B decay backgrounds 
in the B° — > 4>K\ sample are also taken into account. The peaking background fraction 
in the B° — > f (980)Kg sample depends on the PDF used in the fit to the n + 7r~ invariant 
mass distribution, which ignores possible interference between resonant and non-resonant 
amplitudes. We perform a fit to the 7r + 7r~ distribution of a MC sample generated with 
interfering amplitudes and phases for B — > Kixix decays measured from data [19]. The 
observed difference in the signal yield from the true value is taken into account in the 
systematic error determination. We also repeat the fit to the At distribution ignoring the 
contribution of the peaking background. The differences in S and A from our nominal 
results are included in the systematic error. 
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Finally, we investigate the effects of interference between CKM-favored and CKM- 
suppressed B — > D transitions in the /tag final state [28]. A small correction to the PDF 
for the signal distribution arises from the interference. We estimate the size of the correc- 
tion using the B° — > D*~l + v sample. We then generate MC pseudo-experiments and make 
an ensemble test to obtain systematic biases in S and A. We find that the effect on S is 
negligibly small, while a possible shift in A is sizable. 

Various crosschecks of the measurement are performed. We reconstruct charged B meson 
decays that are the counterparts of the B° — > f C p decays and apply the same fit procedure. 
All results for the S term are consistent with no CP asymmetry, as expected. Lifetime 
measurements are also performed for the fcp modes and the corresponding charged B 
decay modes. The fits yield r^o and tb+ values consistent with the world-average values. 
MC pseudo-experiments are generated for each decay mode to perform ensemble tests. We 
find that the statistical errors obtained in our measurements are all consistent with the 
expectations from the ensemble tests. 

For the B° -> <pK° decay, a fit to DS-I alone yields S = -0.68 ± 0.46(stat) and A = 
-0.02 ± 0.28(stat), while a fit to DS-II alone yields S = +0.78 ± 0.45(stat) and A = 
+0.17±0.33(stat). Note that the results for DS-I differ from our previously published results 
S = -0.96±0.50lg;S and A = -0.15±0.29±0.07 [4], as decays B° -> </)K L and <j>K% (K° s -> 
7r°7r°) are included in this analysis. From MC pseudo-experiments, the probability that 
the difference between S values in DS-I and DS-II is larger than the observed difference 
(1.46) is estimated to be 4.5%. A sin 20! measurement with DS-II is performed using 
B° — > J/^Kg (Kg — > 7r + 7r~ or 7T 7r°) and B° — > J/ipK® decays as a crosscheck. Applying 
the same procedure to both DS-I and DS-II, we obtain Sj/^k° — +0.696 ± 0.061 (stat) and 
Aj/^ K o = +0.011 ± 0.043(stat) for DS-I, and Sj/^ K o = +0.629 ± 0.069(stat) and A J/ijjK o = 
+0.035 ± 0.044(stat) for DS-II. The results are in good agreement with each other, and are 
also consistent with SM expectations. As all the other checks mentioned above also yield 
results consistent with expectations, we conclude that the difference in S^k between the 
two datasets is due to a statistical fluctuation. 

A fit to the B° — > fT*°7 (K*° — > K s tt°) sample with both A and S as free parameters 
yields S = — 0.79lo;fo(stat), which is identical with the result of the one parameter fit, and 
A = —0.00 ± 0.38(stat), which is consistent with zero. 

As discussed in Section I, to a good approximation, the SM predicts S = —£f sin 20! for 
the B° -> (j)K°, K + K~K° S , f (980)K° s , r]'K° s , uK° s and K° s ti decays. Figure 3 summarizes 
the sin 20i determination based on our S measurements for these decays. For each mode, 
the first error shown in the figure is statistical and the second error is systematic. For the 
B° — > K + K~Kg decay, an additional systematic error that arises from the uncertainty of the 
CP-even component fraction (lo!oo) i s added in quadrature. We obtain sin 20i = +0.43^o;ii 
as a weighted average, where the error includes both statistical and systematic errors. The 
result differs from the SM expectation by 2.4 standard deviations. 

IV. SUMMARY 

We have performed improved measurements of CP-violation parameters for B° — > (f)K° 
(including both 4>Kg and (pK^), K + K~K$ and 7]' Kg decays, and new measurements for 
B° — > / (980)fr!j, u)K s and K^n decays. These charmless decays are dominated by b — > s 
flavor- changing neutral currents and are sensitive to possible new CP- violating phases. We 
have also measured the time-dependent CP asymmetry in the decay K* ^ (K*° — > -ft^vr ), 
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which is also sensitive to physics beyond the SM. The results for each individual decay 
mode are consistent with the SM expectation within two standard deviations except for the 
B° -> f (980)K° s decay. The combined result for the B° -> (f)K , K + K'K° S , f (980)K° s , 
rj'Kg, ujKg and Kgir decays differs from the SM expectation by 2.4 standard deviations. 
Measurements with a much larger data sample are required to conclusively establish the 
existence of a new CP-violating phase beyond the SM. 
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TABLE I: The event fractions wrong-tag fractions Wi, wrong-tag fraction differences Awi, and 
average effective tagging efficiencies e 1 ^ = e;(l — 2wi) 2 for each r interval for the DS-II. The errors 
for wi and Awi include both statistical and systematic uncertainties. The event fractions are 
obtained from J/ipK^ data. 



1 


r interval 


e/ 


Wl 


Awi 


e eff 


1 


0.000 


- 0.250 


0.397 ±0.015 


0.464 ± 0.007 


±0.010 ± 0.007 


0.002 ± 0.001 


2 


0.250 


- 0.500 


0.146 ±0.009 


0.321 ± 0.008 


-0.022 ±0.010 


0.019 ±0.002 


3 


0.500 


- 0.625 


0.108 ±0.008 


0.224 ±0.011 


±0.031 ±0.011 


0.033 ± 0.004 


4 


0.625 


- 0.750 


0.107 ±0.008 


0.157 ±0.010 


±0.002 ±0.011 


0.051 ± 0.005 


5 


0.750 


- 0.875 


0.098 ± 0.007 


0.109 ±0.009 


-0.028 ±0.011 


0.060 ± 0.005 


6 


0.875 


- 1.000 


0.144 ±0.009 


0.016 ±0.005 


±0.007 ± 0.007 


0.135 ±0.009 



TABLE II: The estimated signal purity and the signal yield iV sig in the signal region for each fcp 
mode that is used to measure CP asymmetries. The result for the B° — > K^ir decay is obtained 
with the sample after flavor tagging but before vertex reconstruction as events that do not have 
vertex information are also used to extract the direct CP violation parameter A. The results for 
the other decays are obtained after flavor tagging and vertex reconstruction. 



Mode 







purity 


iV sig 






-1 


0.63 


139 ± 14 






±1 


0.17 


36 ± 15 


K + K~K° S 




±1 


0.56 


399 ± 28 


fo(980)K° s 




±1 


0.53 


94 ± 14 


V'K° S 




-1 


0.61 


512 ±27 


ujK° s 




-1 


0.56 


31 ±7 




(high-ft s/b ) 


-1 


0.55 


168 ± 16 




(low-ft s/b ) 


-1 


0.17 


83 ± 18 






0.65 


57 ±9 



TABLE III: Results of the fits to the At distributions. The first error is statistical and the second 
error is systematic. We combine B° — > 4>Kg and B° — > 4>K\ decays to obtain S^o and A^o. 



Mode SM expectation for S S A 

lpK° ±sin20i ±0.06 ±0.33 ±0.09 ±0.08 ± 0.22 ± 0.09 

K+K~K° S -sin 20i -0.49 ± 0.18 ± 0.04 -0.08 ± 0.12 ± 0.07 

/ (980)Kg -sin 20i ±0.47 ± 0.41 ± 0.08 -0.39 ± 0.27 ± 0.08 

r]'K° s ±sin20i ±0.65 ± 0.18 ± 0.04 -0.19 ± 0.11 ± 0.05 

uK° s ±sin20i ±0.75 ± 0.64±g;J| ±0.26 ± 0.48 ± 0.15 

K%tt° ±sin20i ±0.30 ±0.59 ±0.11 -0.12 ± 0.20 ± 0.07 

K*° 7 (K*° -> K%ir°) -2(m s /m b )sm2(j) 1 -0.79^ ± 0.10 
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TABLE IV: Summary of the systematic errors on S. 







K+K 




/ (980)^o 


V'K° S 








Vertex reconstruction 


0.01 




0.01 


0.02 


0.01 


0.01 


0.06 


0.02 


Flavor tagging 


0.01 


< 


0.01 


0.01 


0.01 


0.04 


0.02 


0.01 


Resolution function 


0.04 




0.03 


0.03 


0.03 


0.07 


0.05 


0.05 


Physics parameter 


< 0.01 


< 


0.01 


0.01 


< 0.01 


0.01 


0.01 


0.02 


Possible fit bias 


0.01 




0.01 


0.03 


0.01 


+0.01 
-0.10 


0.03 


0.03 


Background fraction 


+0.08 
-0.06 




0.02 


0.05 


0.02 


0.10 


0.02 


0.07 


Background At shape 


0.01 


< 


0.01 


0.04 


< 0.01 


0.02 


0.03 


0.04 


Tag-side interference 


< 0.01 


< 


0.01 


< 0.01 


< 0.01 


0.01 


< 0.01 


< 0.01 


Total 


0.09 




0.04 


0.08 


0.04 


+0.13 
-0.16 


0.10 


0.11 



TABLE V: Summary of the systematic errors on A. 







K + K 


-K° 


/o(980)K° 


V'K° S 






Vertex reconstruction 


0.03 




0.04 


0.04 


0.03 


0.04 


0.04 


Flavor tagging 


< 0.01 


< 


0.01 


0.01 


0.01 


< 0.01 


0.01 


Resolution function 


0.02 




0.02 


0.02 


0.01 


0.04 


< 0.01 


Physics parameter 


< 0.01 


< 


0.01 


< 0.01 


< 0.01 


< 0.01 


< 0.01 


Possible fit bias 


0.01 




0.01 


0.03 


0.01 


+0.01 
-0.03 


0.01 


Background fraction 


0.04 




0.01 


0.06 


0.02 


0.14 


0.02 


Background At shape 


0.03 


< 


0.01 


0.01 


< 0.01 


0.03 


0.01 


Tag-side interference 


0.06 




0.06 


0.03 


0.03 


0.03 


0.05 


Total 


0.09 




0.07 


0.08 


0.05 


0.15 


0.07 
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FIG. 1: The M bc distributions for (a) B° -» <pK%, (c) B° -> K + K~K®, (d) B° -> / (980)lTg, 
(e) B° - (f) B° -> (g) #0 _ (high-7e s/b ), (h) 5° -> K°vr (low-^ s/b ), and 

(i) B° -> ET*°7 within the A_E signal region and (b) the p^ 118 distribution for B° — > c/>-KT£. Solid 
curves show the fit to signal plus background distributions, and dashed curves show the background 
contributions. 
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FIG. 2: The asymmetry, A, in each At bin with < r < 0.5 (top) and with 0.5 < r < 1.0 (bottom) 
for (a) B° -» 0K°, (b) B° -» K + K~K° S , (c) 5° -► / (980)lifg, (d) B° -► r/.Kg, (e) 5° -> wlfg, 
(f) i? — > -f^gvr (for the high-7£ s /b region only), and (g) B° — > K*°^f. The solid curves show the 
result of the unbinned maximum- likelihood fit. The dashed curves show the SM expectation with 
sin 20i = +0.73 (5 = for B° -» K*°>y) and A = 0. 
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Belle Results on Time-dependent CP Violation in b— »s 
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FIG. 3: Summary of sin2</>i measurements performed with B decay governed by the b — > sqq 
transition. The world-average sin2c/>i value obtained from B° — > J/t/jK° and other related decay 
modes governed by the b — > ccs transition [6] is also shown as the SM reference. 
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